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ABSTRACT 

We report the results of an X-ray proper motion measurement for the NW 
rim of SN1006, carried out by comparing Chandra observations from 2001 and 
2012. The NW limb has predominantly thermal X-ray emission, and it is the only 
location in SN 1006 with significant optical emission: a thin, Balmer-dominated 
filament. For most of the NW rim, the proper motion is ~ 0.30" yr~^, essentially 
the same as has been measured from the Ha filament. Isolated regions of the NW 
limb are dominated by nonthermal emission, and here the proper motion is much 
higher, 0.49" yr~^, close to the value measured in X-rays along the much brighter 
NE limb, where the X-rays are overwhelmingly nonthermal. At the 2.2 kpc 
distance to SN1006, the proper motions imply shock velocities of ~ 3000 kms~^ 
and ~ 5000 kms~^ in the thermal and nonthermal regions, respectively. A lower 
velocity behind the Ha filament is consistent with the picture that SN 1006 is 
encountering denser gas in the NW, as is also suggested by its overall morphology. 
In the thermally-dominated portion of the X-ray shell, we also see an offset in the 
radial profiles at different energies; the 0.5-0.6 keV peak dominated by O VII is 
closer to the shock front than that of the 0.8-3 keV emission — due to the longer 
times for heavier elements to reach ionization states where they produce strong 
X-ray emission. 
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Introduction 



The supernova of SN 1006 was observed widely in Asi a and Europe, an d appears to 
have been the brightest naked-eye supernova ever recorded ( IStephensonll2010l ). Its remnant 
appears in catalog listings as SNR G327.6+14.6, though in this paper we adopt the common 
parlance of using "SN1006" to refer to the remnant. SN 1006 stands out ar nong historical 
remn ants as largest (30'= 19 pc diameter), one of the closest, at 2.2 kpc (j Winkler et al. 
20031 ). farthest above the Galactic plane (550 pc), and least obscured {Nh ~ 5 x lO^'' cm~^, 
all of which make it amenable to study in X-rays and other bands. However, it is also the 
faintest of the historical remnants, in total flux and especially in surface brightness, so deep 
observations are required for detailed study. Since SN 1006 is both a prototypical Type la 
SNR, and since it shows some of the best and cleanest evidence for electron acceleration to 
TeV energies in shocks, its close examination can further our understanding of both SNR 
dynamics and particle acceleration in general. 

The overall structure of SN 1006 is that of a nearly circular limb-brightened shell with 
strong bilateral symmetry; its N E and SW limbs are far brighter, in both radio and X-rays 
than other portions of the rir n (ICassam-Cheriai et al.l l2008l : iDyer et al.l 120091 : iMiceli et al. 
20091 . and references therein). iKoyama et al.l (119951 ) found that the X-rays from the bright 
NE and SW limbs are not only harder than those from the rest of SN 1006, but also that 
their spectrum is a featureless power-law, a result that provided the first clear evidence 
for diffusive acceleration of charged particles and cer nented the long- suspected link between 
supernova remnant (SNR) shocks and cosmic rays (iReynoldd Il996l ). Furthermore, the X- 
ray morphology along the NE and SW rims matches the radio images in exquisite detail 
(jWinkler fc Long! 119971 ) , further confirming a common synchrotron origin for emission in 
both widely separated bands. 

The NW and SE limbs of SN 1006 are dramatically different; the shell in both the NW 
and SE is much fainter and less distinct in the X-ray and radio bands, and the X-ray spec- 



^Based on observations made with NASA's Chandra X-ray Observatory, NASA's Chandra Observatory is 
operated by Smithsonian Astrophysical Observatory under contract # NAS83060 and the data were obtained 
through program GOl-12115. 
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trum is soft and thermal, dom inated by lines from He-like and H-like oxygen between 0.5 



and 0.8 keV (ILong et al.ll2003[ ). The NW limb is unique in that only here can significant 



Schweizer &; Laskei 


19781: 


supernova ( 


Kirshner & Chevalier 



19781 ). Balmer-dominated, or "nonradiative" filaments arise when a fast shock expands into 
a low-density, partially neutral environment. The Balmer emission has two components: 
neutral H atoms that enter the shock can be excited, and subsequently decay to produce 
narrow lines (with a width characteristic of the pre-shock temperature); or they can undergo 
charge exchange with hot protons to pr oduce broad lines whose width is closely related to 
the p ost-shock proton temperature (e.g.. lChevalier et al.lll980l : iGhavamian et al.ll2002l : iHeng 



2010l ). Since the lifetime of neutral atoms in the hot post-shock environment is very short, 
the Balmer filaments can occur only immediately behind the shock, and thus delineate the 
current position of the shock. 



Hesser & van den Berg] 


1 ( 


1981) 


recently by 


Long et al. 


(1988|) c 



CCD images taken over an 11-year baseline, obtained M Hr> = 0.280 ± 0.008" yr \ Based on 



the Balmer profiles and models for the shock conditions, IGhavamian et al.l (120021 ) obtained a 
value of Vs = 2890 ± 100 kms~^ for the shock velocity in the NW, and combining this value 
with the proper-motion measurement yields a geometrically determined distance of 2.2 ±0.1 
kpc to SN 1006. Having a relatively precise distance of course means that proper motions 
measured by any means, anywhere around the SN 1006 shell, translate directly into a shock 
velocity at that point, which is si gnificant becau se Vg is the most crucial parameter in the 



theory of shock acceleration (e.g.. lReynoldsll2008l ). 



A useful way of expressing the proper motion is through the expansion parameter m: 
the power- law index in i? oc t™, where t is the age of the remnant. This parameter can be 
interpreted as the ratio of current expansion rate divided by the mean rate over the remnant's 
lifetime, m = lit/ 9, where 9 is the angular radius. The expansion index gives some insight 
into a remnant's evolutionary state, e.g., m = 1 for free expansion, or m = 0.4 for Sedov 
expansion. 

Proper-motion measurements of SN 1006 have been made in various bands with dif- 
ferent precision. Radio measurements of the global expansion gave mradio = 0.48 ± 0.13 



^In addition to the NW filam ents, there is exceedingly faint and more diff use Balmer emission surrounding 
almost all of the SN 1006 sheU (jwinkler fc Lon j|l997t [winkler et al.ll2003[ ). 
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( iMoffett et al.lll993l ). c onsiderably larger t han the value of mna = 0.34 ± 0.01 for the NW 
optical filament alone ( IWinkler et al.ll2003[ ). Since the stronger optical emission in the NW 
indicates a higher external density, it is not surprising that the expansion should be slower 
than the average around the shell, since the higher the density, the more rapidly the expan- 
sion will decelerate. (If one assumes pressure equilibrium, which may pertain in SN 1006 
but is not assured, the preshock density will scale inversely as the square root of the shock 
velocity.) In a recent paper, we obtained an X-ray measurement of the proper motions along 
the E-NE limb of SN 1006, using a pair of Chandra ACIS-S observations from 2000 and 
2008, and found u = 0.48 ± 0.04" yr-\ or m = 0.54 ± 0.05, and Vs = 5000 ± 400 kms"^ 



flKatsuda et al.ll2009[ ) 



In this paper, we present the first X-ray measurement of the proper motion along the 
NW limb of SN1006, using a pair of deep Chandra ACIS-S observations from 2001 and 
2012. These new measurements cover much the same region as the optical proper motion 
measurement by I Winkler et al.l (120031 ) . and the results are consistent. In addition to the 
thermal X-ray emission that dominates most of the NW limb, there are two small regions 
with a much harder X-ray spectrum, free of emission lines and consistent with nonthermal 
synchrotron emission. We have also measured t he motions in these regions and find these 
to be consistent with our earlier measurements (iKatsuda et al.l 120091 ) for the much brighter 
synchrotron-dominated NE rim. The 2012 April observation of the NW limb represents the 
first in a series of deep ACIS images covering the entirety SN 1006, being carried out as a 
Chandra Large Project. Further results will be reported in subsequent publications. 



2. Observations and Results 

The initial Chandra observation of the NW limb was carried out with the ACIS-S in 



2001 April and reported by iLong et al.l (120031 ). In order to measure proper motions in the 



NW, we carried out a repeat observation in 2012 April; the roll angle and exposure time were 
almost identical for both observations, and the positions were similar, with the second-epoch 
one displaced slightly to the NW in order to include more of the pre-shock background. The 
two observations are summarized in Table [H In both cases we reprocessed the level-1 event 
files with CIAO ver. 4.4 and CALDB ver.4.5.1. 

A three-color image of a portion of the NW limb of SN 1006 as seen with Chandra in 
2012 is shown in Figured], where red, green, and blue correspond to 0.5-1 keV, 1-2 keV, and 
2-7 keV bands, respectively. Most of the X-ray emission in the NW is thermal, and appears 
reddish in this image, but there are also a few fainter nonthermal, synchrotron-dominated 
regions that appear more bluish. The difference in spectral character between the thermal 
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and nonthermal regions is most apparent in spectra from these regions, shown in Fig. [2l 
For the sharp, bright shock front in the NE, an outward displacement was evident in 



a sim ple difference image between the images at epochs eight years apart (iKatsuda et al. 



20091 ). In the NW, however, the surface brightness is much lower and the shock front less 
sharp, so the expansion is less obvious. In order to measure it, we have integrated along the 
shock front to give radial profiles in three thermal regions where the shock is well defined, 
labeled Tl, T2, and T3, and also in two nonthermal regions, labeled NTl and NT2. 

As described in Sections 2.2 and 2.3, we have measured the proper motions in different 
regions in different energy bands. In order to produce fiux-calibrated images, it was necessary 
to consider the variation in effective exposure time with both energy and position on the 
detector. We created exposure maps at monochromatic energies: 0.6 keV, 0.55 keV, 1.2 keV, 
and 2keV for 0.5-1 keV, 0.5-0.6 keV, 0.8-3 keV, and 1-8 keV bands, respectively — energies 
that roughly correspond to mean photon energies in the energy bands of interest, given the 
spectral shape in each region (Fig. |2]). 



2.1. Image Registration 

Before measuring the proper motion, we first registered the images using as fiducials four 
bright point sources, all with significance levels greater than 10 a and point-spread-function 
(PSF) sizes less than 2" in both imagesjf] We used the CIAO tool wavdetect to obtain 
positions for the sources at the two epochs, which we summarize in Table [21 the sources are 
also indicated in the cyan circles on Fig. [T] Sources PI and P3 correspond to faint [R ~ 19) 



star-like objects in the NOMAD catalog, (jZacharias et al.ll2005l ). both with negligible proper 
motion, while the P2 and P4 sources seem to have no optical companion0 All four sources 
have a hard spectrum, and all are likely to be background AGNs; we assume that all have 
negligible proper motions. The offset between the two images, obtained from a simple mean 
of the position differences listed in the two right-hand columns of Table [2], are 0.06" and 
0.18", in R.A. and Decl., respectively. Such offsets are consistent with the stated absolute 
astrometric accuracy for Chandra and the ACIS-S arraylfl After correcting for these offsets. 



•^Here PSF size is defined as a radius in which 39% of total source photons are enclosed. 

''^Source P2 is located 2.5" NW of another faint object in the NOMAD catalog, but since the position dif- 
ference is well outside the error ellipses from both ACIS observations, this is unlikely to be a true association. 
In any case, this object also has no proper motion in the NOMAD catalog. 

^http://cxc.harvard.edu/cal/ASPECT/celmon/. 
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the RMS residuals in position difference between the two epochs for the four point sources 
are 0.37" in R.A. and 0.27" in Decl. In our proper motion analysis that follows, we have 
apphed the mean offsets, and have included both statistical uncertainty and uncertainty 

in the registration (for the radial direction), calculated as ^(0.37" sin 6^)^ + (0.27" cos 6')^, 
where 9 is the azimuth angle for each area, measured counterclockwise from North. 



2.2. Thermal Regions 



We selected regions Tl, T2, and T3 both because they coincide with regions studied 
previously and because the shock front is sharply defined there. Region Tl includes a small, 
bright X-ray knot, located within an Ha bubble a h ead o f most of the NW shock, and also 
includes the bright X-ray knot where IVink et al.l ( 120031 ) directly measured both electron 
and ion temperatures using XMM-Newton reflection grating spectra. They demonstrated 
that Tg remains much lower than Tj well behind the shock, and also that non-equilibrium 
ionization conditions prevail. Regions T2 and T3 trace the shock front delineati ng the NW 



limb where proper motion s of the Ho; filaments have been precisely measured (ILong et al. 
19881 : Iwinkler et al.ll2003l ). 



In Figure [31 we plot the projected one-dimensional profiles for epochs 2001 and 2012, 
which clearly show outward motion in all three regions. The profiles have been generated 
from vignetting-corrected images in the 0.5-1 keV band, which contains a large fraction of the 
total X-ray emission. To quantitativel y measure the shifts between epochs, we have followed 
the approach we used previously (e.g.. lWinkler et al.ll2003l : iKatsuda et al.ll2009l ). Briefly, the 
method is to minimize the value for the difference between the second-epoch profile and 
the shifted first-epoch one, as a function of the amount of the shift. When calculating the 
values, we use data points within the vertical dashed lines indicated in each panel of Figure [3l 
so that we can concentrate on features of interest. The results of these measurements for all 
three regions are given in Table [31 

The X-ray proper motion we measure in region T2 is consistent, w ithin the uncertainty , 
with the optical measurement for the same region (0.280"±0.008"yr~^: IWinkler et al.ll2003l ). 



The proper motions in regions Tl and T3 are slightly higher than in T2, b ut all are con 



sistent within the uncertainties. Taking the distance to SN 1006 as 2.2 kpc (IWinkler et al. 



20031 ) ■ the expansion velocity for Tl is 3300 ± 200 ± 300 kms~^ (statistical and registra- 
tion uncertainties, respe ctively), in r e asona ble agreement with the velocity inferred from the 
X-ray line width, where IVink et al.l ( 120031 ) found that the velocity is > 4000 kms~^, if no 
significant temperature equilibration has taken place, but could be as low as ~ 3000 kms~^ 
with some adiabatic cooling. 
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The first radio measurement of tlie SN 1006 expansion was by lMoffett et al.l ( 119931 ). who 
measured a global average expansion rate of 0.44"±0.13" yr^^ over a baseline of only 8 years. 
They also give local rates in four broad azimuthal sectors, but the NW rim was too indistinct 
to give any measurement at all there. However, iMoffett et al.l ( 120041 ) have since obtained a 
third-epoch radio image of SN 1006, and with a longer baseline obtained expansion rates in 
the NW that are consistent with the optical values. 

In Table [31 we also give the local values for the expansion index, m. Calculating m = 
fj,t/9 requires an angular radius, 9, and to define this we have taken the expa nsion center to 
be R.A . (J2000.) = 15:02:54.9, Decl. (J2000.) = -41:56:08.9, as determined bv lXatsuda et al. 
( 120091 ) from the i^O^^THRI X-ray image mosaic of the entire remnant. 



For Area T2, we have also examined radial profiles and proper motions in distinct energy 
bands. Of the three thermal regions, T2 is the best suited for such a study because it has 
the sharpest shock front and the best photon statistics. As shown in Figure H] (left), the 
radial profiles vary with X-ray energy; the peak of the 0.5-0.6 keV (a band dominated by 
O VII triplets) profile is located just behind the Ha filament, whereas the 0.8-3 keV peak 
is broader and shifted by ~ 15" further behind the shock. This must be primarily due to 
the evolution of the ionization to produce He-like species of elements h eavier than o xygen , 
as shown in the profiles and spectra along the NW rim that appear in iLong et al.l ( 120031 ). 
Ionization of higher-Z metals, as well as the gradual increase in the electron temperature 
behind the shock, both make the X-ray spectrum harder toward the center of the SNR. We 
have measured proper motions in these two energy bands, using the one- dimensional profiles 
shown in Figure H] (center and right), and find no significant energy dependence; the values 
are given in Table [31 



2.3. Nonthermal Regions 

We have also carried out proper-motion measurements for the small regions in the NW 
that are dominated by nonthermal X-rays, i.e., areas NTl and NT2 in Figure [H Since the 
X-ray emission in these areas is harder, we have used a higher and broader energy band, 
1-8 keV, to measure the proper motion. The radial profiles in Figure [51 show clear outward 
motions that are significantly larger than those found for the thermal regions. Measurements 
using the same minimization technique yield almost identical values of 0.48 to 0.49" yr~^ 
for both NTl and NT2, both ~ 50% higher than the values obtained in the thermal regions. 
Furthermore, the proper motions in NW regions NTl and NT2 are consistent with our 
previous measurements of 0.48 ± 0.04" yr~^ measured at more than a dozen locations along 
the shock in the much brighter synchrotron-dominated NE region of SN 1006 (iKatsuda et al. 



- 8 - 



2009). 



Discussion 



We have confirmed that for the thermal-dominated NW shock, the proper motions 
measured in X-rays are consistent with those measured optically for the Balmer-dominated 
filaments that delineate the shock. While not unexpected, this is the first time that proper 
motions have been measured in both X-ray and optical bands for the same region in SN 1006. 
This is noteworthy because for other young remnants, past measurements in different bands 



have sometimes produced widely discrepant results (e.g., see the summary in iKatsuda et al. 



20091 ). though these may well have been due to comparing local measurements with global 
average s. For the Tyc h o SN R, at least, recent high-resolution measurements with Chan- 
dra bv iKatsuda et al.l (120101 ) are in reasonable agreement with radio measurements by 
Reynoso et al. ( 19971 ) around the entire shell. We have now seen that for SN 1006, X-ray and 
optical measurements along the NW rim are in excellent agreement. But as the capability for 
detailed measurements becomes available, it is important to keep in mind that the features 
themselves can evolve, and that new features can appear while others disappear, giving the 
illusion of proper motions in snapshot images taken many years apart. 

Our confirmation that the shock velocity in the thermal-dominated NW limb is slower 
than at other locations around the SNR shell provides direct evidence that the SNR shock 
is interacting with a dens er ambient medium in the NW region, as has previously been sug- 
gested by several studies. iHeng et al.l (120071 ) found that the Balmer-dominated Ha filaments 
in the NW ( which have significan tly higher surface brightness than optical emission elsewhere 
in SN1006, IWinkler et al.ll2003l ) are consistent with a pre-shock density of 0.15-0.3 cm~'^, 
higher by a factor of several from that around the rest of the remnant. This density is simi- 
l ar to that inferred from X-ray spectr a from Chandra ( Long et al.l 2003 ) and XMM-Newton 



(Acero et al.l 120071: iMiceh et al 



20091). as well as from FUSE UV spectra (IKorreck et al. 



2004 ). Most recently, IWinkler et al.l (120121 ) report the detection of 24/im infrared emission 
arising from warm dust grains in SN 1006. This emission is seen only along the NW rim, 
and Spitzer spectra of the warm dust is consistent with a post-shock density of 1 cm~'^. 

The shock in the NW of SN 1006 is apparently encountering a localized region of higher 
density than that surrounding the rest of the remnant, possibly representing t he start of an 



encounter with the extended H I concentration just beyond the rim found by iDubner et al. 



(I2OO2I ). As expected where the pre-shock density is higher than average, the rim in the NW is 
indeed flatter than elsewhere, with a radius up to 15% smaller there than the average around 
SN 1006. Since the NW filament lies at a radius only slightly smaller than the remnant mean. 
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this slower expansion cannot have been going on for very long, and so the shock must have 
encountered the denser region fairly recently. 

Another point of view is to consider the expansion index m. In the E and NE, m ~ 0.54, 
indicating that the expansion is well short of reaching the Sedov pha se (m = 0.4) in thes e 



regions, consistent with expansion into a very low-density environment (IKatsuda et al.ll2009l ). 
But in the NW, where the radius is only slightly smaller than elsewhere, m ^ 0.38, consistent 
with Sedov expansion, and thus with a higher density. 

The small sections of the NW limb where nonthermal X-ray emission dominates have 
proper motions, and hence shock velocities, significantly higher than those in the thermal 
regions. As noted above, the measured values of ~ 0.5" yr~^ for the nonthermal regions 
NTl and NT2 give shock velocities ~ 5000 kms~^ at a dista nce of 2.2 kpc, simila r to those 
all along the nonthermal, synchrotron-dominated NE limb (IKatsuda et al.ll2009l ). Thus it 
appears likely that in SN 1006, nonthermal-dominated regions are commonly associated with 
fast shocks of ~ 5000 kms~^. Higher shock speeds than in the thermal regions indicate a 
lower pre-shock density for the the nonthermal regions, suggesting the existence of low- 
density pockets within the generally higher ISM density to the NW, or other small scale 
inhomogeneities. 

Higher shock speeds in the nonthermal regions are also consistent with the theoretical 
view that faster shocks can enhance synchrotron X-ray emission by boosting the roll-off 
frequency z/roUoff- Diffusive shock acceleration theory predicts a power-law electron energy 
distribution with an exponential cutoff at some charac teristic energy, and ^^miinff is the peak 
frequency emitted by electrons with that energy (e.g., lEUison et al.ll2000l ). iReynoldsl (120081 ) 
has shown that froiioff oc when the maximum energy of accelerated particles is limited by 
the SNR age, or z/j-oiioff oc when synchrotron losses limit the maximum energy. 

Still, a shock speed of ~ 3000 kms~^ should be high enough for substantial particle ac- 
celeration to X-ray-emitting energies. The strength and orientation of the upstream magnetic 
field, as well as the electron diffusion coefficient, also play important roles in determining 
acceleration efficiency, and could inhibit particle acceleration in the NW while enhancing it 
in the NE and SW. In addition, a partially neutral pre-shock medium, a condition that the 
existence of bright Balmer filaments indicate must pertain along most of the NW limb, may 
inhibit shock acceleration (IRaymond et al.ll201ll : iBlasi et al.ll2012l ). In particular, Blasi et al. 
find that for shock speeds below 3000 kms~^, the energy spectrum for accelerated particles 
can steepen markedly at low energies, reducing their X-ray synchrotron emission compared 
to that from faster shocks. 



Finally, we note that the SE portion of the remnant has predominantly thermal X- 
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ray emission (iRothenflug et al.ll2004l ). but resembles the nonthermal regions in terms of the 
observed local remnant radius (which implies the mean shock velocity over the lifetime of 
the remnant) and in t he fai ntness of its Ha emission (which implies the ambient density). 



Cassam-Chenai' et al.l ( 120081 ) have attributed the dramatic difference largely to a far lower 
efficiency for particle acceleration to the SE than in the NE or SW. Measuring proper 
motions along the SE rim should be an important diagnostic for understanding the nature 
of particle acceleration in SN 1006. 



4. Summary 



We have measured X-ray proper motions in the NW limb of SN 1006, using the first 
of our new Chandra Large Project observations. For the thermally-dominated X-ray re- 
gions along most of the NW limb, we find proper motions of 0.28-0.35" yr~^, consistent 
(wi thin the uncer t aintie s) with optical measurements for the Balmer filaments, 0.28" yr~^, 
by IWinkler et al.l (120031 ). Even along the thermal-dominated NW rim, however, we find 
two small, isolated regions with harder, non-thermal X-ray emission. Proper motion mea- 
surements fo r these regions are mu ch higher: /ix ~ 0.49" yr~^, almost identical to the values 
measured by lKatsuda et al.l (120091 ) along the synchrotron-dominated E to NE rim of SN 1006. 
We attribute the slower motion along most of the NW rim to the relatively recent encounter 
of the SNR shock with a higher density pre-shock sheet. However, the existence of small re- 
gions with non-thermal X-ray spectra and much higher shock speeds in close proximity to the 
more extensive regions with thermal spectra and slower shock speeds suggests low-density 
pockets within this denser sheet, or other small-scale inhomogeneities. 
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Table 1. Chandra ACIS-S Observations of SN 1006, NW 



ObsID 


RA(J2000.) 


Dec(J2000.) 


Roll 


Obs. Date 


Exposure 


PI 


1959 


15:02:22.5 


-41:47:08 


30.2° 


2001 Apr 26 


89.0 


Long 


13737 


15:02:15.9 


-41:46:10 


31.7° 


2012 Apr 20 


87.1 


Winkler 



Table 2. Positions of Fiducial X-ray Point Sources 



Source 




2001 






2012 




Difference'' 


R.A. 


Decl. 


PSF size* 


R.A. 


Decl. 


PSF size* 


R.A. 


Decl. 


PI 


15:02:07.010(0'.'13) 


-41:51:26.70(0'.'15) 


I'.'S 


15:02:07.022(0'.'12) 


-41:51:26.13(0'.'15) 


1'.'9 


-or.' IS 


0'.'57 


P2 


15:02:30.468(0'.'26) 


-41:42:35.20(0'.'36) 


1'.'3 


15:02:30.470(0'.'20) 


-41:42:35. 12(0'.'21) 


1'.'2 


-0'.'02 


0'.'08 


P3 


15:02:34.514(0'/22) 


-41:42:00.53(0'.'18) 


l'/8 


15:02:34.54(0'.'15) 


-41:42:00.69(0'.'09) 


1'.'7 


-0'.'29 


-o'.'ie 


P4 


15:02:40.555(0^'13) 


-41:45:09.02(0'.'12) 


0'/8 


15:02:40.495(0'.'17) 


-41:45:08.79(0'.'13) 


1'.'3 


0'.'67 


0'.'23 



^Radius in which 39% of total source photons are enclosed. 
Positive values indicate apparent motions to the west and north in R.A. and Decl., respectively. 



Note. — 



Values in parentheses are 1-a statistical uncertainties. 
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Table 3. Summary of Proper-Motion Measurements 



Region 


Energy Band 
(keV) 


Radius 
(') 


Proper Motion** 
("yr-i) 


Expansion Index** 


Velocity**'^ 
( kms~^) 


(Xi/)min 


Tl 


0.5-1 


14.1 


0.32±0.03±0.03 


0.38±G.03±G.04 


3300±200±3GG 


0.99 


T2 


0.5-1 


13.3 


0.28±0.02±0.03 


0.34±G.03±G.G4 


2800±200±3GG 


0.58 


T2b 


0.5-0.6 


13.3 


0.28±0.03±0.03 


0.35±0.04±0.04 


2900±300±300 


0.76 


T2c 


0.8-3 


13.3 


0.28±G.06±G.03 


0.35±G.G8±G.G4 


2900±600±3GG 


0.94 


T3 


0.5-1 


14 


0.35±G.05±G.03 


0.42±G.G6±G.G3 


3700±500±3GG 


0.84 


NTl 


1-8 


14.7 


0.48±G.05±0.03 


0.54±0.G6±G.04 


5000±500±300 


1.32 


NT2 


1 8 


15.2 


().49±().U4±(),U3 


0.53±().U5±().03 


3000±4()()±2()() 


1.13 



"^The first and second error terms represent 1-a statistical and registration uncertainties, respectively. 
^Assumes a distance of 2.2 kpc. 
•^Minimum per degree of freedom. 
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Fig. 1. — Chandra three-color image of the northwestern hmb of SN 1006. Red, green, 
and blue correspond to fluxes in the 0.5-1 keV, 1-2 keV, and 2-7 keV bands, respectively. 
The image is binned by 1" and has been smoothed by a Gaussian kernel of a = 6". The 
intensity is scaled as the square root of the count rate. The fields of view of the Chandra 
observations of the northwestern limb (ObsIDs 1959 and 13737 covering relatively inner and 
outer regions, respectively) are shown as white boxes with ACIS chip identifications, and the 
four fiducial sources (PI - P4 from west to east) used to register the images are indicated by 
the cyan circles . The regions where we extract radial profiles are shown in white (thermal 
dominated) and yellow (nonthermal dominated); note the spectral difference between the 
prominent features in the thermal vs nonthermal regions. The field measures 17' x 13' and 
is oriented north up, east left. 
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Fig. 2. — Chandra ACIS spectra extracted from the thermal and non-thermal regions indi- 
cated in Fig. 1. For clarity, spectra from different regions have been scaled by the factors 
indicated. The spectra have been extracted not the entire the regions shown in Fig. 1, but 
instead from the filamentary (T2, T3, and NT2) or knotty (Tl and NTl) features that we 
used for the proper- motion measurement. 
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Fig. 3. — Radial profiles extracted from the thermally-dominated areas shown in Figure [H 
Data points in black and red represent the 2001 and 2012 epochs, respectively. The intensity 
in 2012 is scaled to equalize that in 2001 by factors of 0.93, 1.01, and 1.11 for Tl, T2, and T3, 
respectively. The dashed vertical lines demarcate the regions used in the measurement of 
the proper motions. The profile labeled Tl (left panel) covers both the Tl and NTl regions 
in Figure [H the smaller peak represents emission from the NTl region, which is expanding 
faster than the thermal-dominated Tl peak. Figure [5] (left) shows the profile from exactly 
this same area, but for higher energy X-rays. 
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Fig. 4. — Left: Radial profiles in the T2 region. Black, red, and blue represent Ha, 0.5- 
0.6 keV, and 0.8-3 keV emission, respectively. The X-ray data taken in the two epochs are 
combined after taking account of proper motions which amounts to 3 arcse c between the 
two o bservations. The Ha profile is shifted by an expected proper motion (jWinkler et al. 
20031 ) at the X-ray observation. Center and right: Same as Figure |3l but for different energy 
bands. The intensity in 2012 is scaled to equalize that in 2001 by factors of 1.01 and 0.93 
for 0.5-0.6 keV and 0.8-3 keV, respectively. 
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Fig. 5. — Same as Figure |3l but for the NTl and NT2 regions, using the 1-8 keV band. The 
intensity of 2012 is scaled to equahze that of 2001 by factors of 1.09 and 0.94 for NTl and 
NT2, respectively. The profile labeled NTl (left panel) is from exactly the same area as the 
one labeled Tl in Figured] above 1 keV, emission from the outer, non-thermal region (NTl) 
is stronger than that from the inner, thermal one (Tl). 



